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ABSTRACT 


The  evolution  of  the  frequency-directional  wave  spectrum,  E(/j9),  across  the  inner 
continental  shelf  and  beach  was  examined  with  measurements  collected  at  the  U.S.  Army  Corps 
of  Engineer’s  Field  Research  Facihty  during  the  recent  SandyDuck  experiment.  Arrays  of  bottom 
pressure  sensors  were  deployed  on  the  shelf  in  20  m  depth  and  on  the  beach  in  depths  ranging 
from  2  -  5  m.  These  arrays  were  complemented  by  a  directional  wave  buoy  in  20  m  depth  and 
an  array  of  pressure  sensors  in  8  m  depth  maintained  by  the  U.S.  Army  Corps  of  Engineers.  A 
preliminary  analysis  of  these  data  is  presented  here  focused  on  four  case  studies  that  illustrate  the 
observed  wave  shoaling  evolution  in  both  non-breaking  and  breaking  conditions.  Estimates  of 
E(/,6)  extracted  from  array  cross-spectra  at  six  cross-shore  locations  are  compared  to  predictions 
of  linear  refraction  theory.  The  present  observations  support  conclusions  from  previous  studies 
that  the  cross-shore  evolution  of  dominant  wave  propagation  direction  is  well  described  by  linear 
refraction  theory.  Observations  of  harmonic  peak  development  at  directions  aligned  with  the 
dominant  waves  are  consistent  with  theoretical  wave-wave  interaction  rules  and  previous 
observations.  In  both  non-breaking  and  breaking  conditions,  the  observed  E(f,9)  are  directionally 
broader  than  predicted.  In  contrast  to  previous  observations  on  a  barred  beach,  the  present 
observations  on  a  planar  beach  do  not  show  a  dramatic  broadening  of  directional  wave  spectra 
in  the  surf  zone. 
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I.  INTRODUCTION 


Accurate  predictions  of  the  nearshore  transformation  of  ocean  waves  are  vital  to 
naval  littoral  operations  (e.g.,  mine  warfare,  amphibious  landings,  swimmer  operations) 
and  civilian  coastal  erosion  studies.  As  waves  propagate  from  the  open  ocean  onto 
shallow  coastal  shelves  and  ultimately  break  on  beaches,  linear  and  nonlinear  effects 
change  their  amphtudes  and  directional  characteristics.  The  present  study  is  focused  on 
the  still  poorly  understood  transformation  of  the  directional  properties  of  waves 
propagating  across  the  irmer  shelf  and  beach. 

Linear  shoaling  of  waves  is  well  understood.  Reduction  of  the  group  speed  in 
shallow  water  generally  causes  an  amplification  of  wave  amplitudes  on  beaches.  Waves 
arriving  at  an  oblique  angle  relative  to  the  shore  are  refracted  toward  normal  incidence 
while  for  large  incidence  angles  their  amplitudes  are  reduced  (e.g..  Kinsman,  1965). 

Close  to  shore,  strong  nonlinearities  and  wave  breaking  affect  the  frequency 
spectra  and  directional  properties  of  waves.  Nonlinear  interaction  between  two  primary 
wave  components  with  frequencies  /j  and  ^  and  vector  wavenumber  ki  and  ki  excite  a 
secondary  wave  component  with  the  sum  frequency  (fi+f2)  and  wavenumber  (ki+ka).  If 
the  incident  wave  spectrum  is  narrow  (i.e.,/j«/2,  ki^kz),  these  secondary  waves  are 
harmonic  components  with  approximately  twice  the  frequency  of  the  incident  waves  and 
propagate  in  approximately  the  same  direction  as  the  incident  waves  (Freilich  et  al.  1990; 
Berbers  et  al.  1998).  Subsequent  interactions  between  the  primary  waves  and  the  newly 
formed  secondary  waves  excite  higher  harmonic  components  (3/i,  4/j,  ...)  all 
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approximately  directionally  aligned  with  the  incident  waves. 

Field  observations  of  the  shoaling  evolution  of  directional  wave  spectra  are  scarce. 
Freilich  et  al.  (1990)  presented  high  resolution  directional  spectrum  estimates  obtained 
from  two  arrays  of  pressure  sensors  located  well  outside  the  surf  zone  on  a  nearly  planar, 
natural  beach.  These  observations  show  the  expected  refraction  toward  normal  incidence 
of  incident  swells  and  energy  transfers  to  directionally  aligned  harmonics.  In  bimodal 
spectra,  with  two  distinct  wave  systems  arriving  at  the  beach  from  different  incidence 
angles,  the  theoretically  expected  excitation  of  waves  at  the  sum  frequency  and  vector 
wavenumber  was  observed.  Berbers  et  al.  (1998)  presented  estimates  of  mean  directions 
and  directional  spreading  parameters  extracted  from  cross-shore  transects  of  co-located 
pressure  sensors  and  current  meters  on  barred  and  non-barred  beaches.  These 
observations  show  that  although  the  mean  wave  direction  does  not  appear  to  be  strongly 
influenced  by  wave  breaking,  the  directional  spreading  of  waves  increases  significantly 
in  the  surf  zone  on  both  barred  and  non-barred  beaches. 

In  the  present  study,  a  preliminary  analysis  is  presented  of  a  more  extensive 
shoaling  wave  data  set  acquired  dming  the  SandyDuck  experiment  near  Duck,  NC.  Six 
high  resolution  arrays  of  bottom  pressure  sensors  were  deployed  for  four  months  along 
a  transect  extending  from  the  inner  shelf  through  the  surf  zone.  These  measurements 
provide  unique  detailed  estimates  of  the  shoaling  and  breaking  transformation  of 
directional  wave  spectra.  The  field  experiment  and  data  analysis  techniques  are  described 
in  chapter  11.  Four  case  studies  spanning  a  range  of  conditions  (i.e.,  nonbreaking  and 
breaking  waves,  narrow  swell  spectra,  and  broad  wind  seas)  are  described  in  chapter  III. 


2 


The  results  are  summarized  in  chapter  IV. 
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n.  FIELD  EXPERIMENT  AND  DATA  ANALYSIS 


Detailed  measurements  of  the  evolution  of  frequency-directional  wave  spectra 
across  the  iimer  continental  shelf  and  beach  were  collected  during  the  SandyDuck 
experiment  at  the  U.S.  Army  Corps  of  Engineer’s  Field  Research  Facility  located  near 
Duck,  NC,  on  a  straight  barrier  island.  This  site  is  fully  exposed  to  the  Atlantic  Ocean 
and  features  an  80  km  wide,  shallow  (20  -  50  m  depths)  shelf.  Arrays  of  pressure  sensors 
were  deployed  on  the  iimer  shelf  in  20  m  depth  (5  km  from  shore)  and  on  the  beach  in 
depths  ranging  from  2  -  5  m  (Figure  1).  The  experiment  was  conducted  in  collaboration 
with  Dr.  Steve  Elgar  (Washington  State  University)  and  Dr.  R.  T.  Guza  (Scripps 
Institution  of  Oceanography).  High-quality  wave  data  were  collected  at  aU  instrument 
locations  during  a  four  month  period  (August  -  November  1997). 

The  20  m  depth  array  consisted  of  nine  pressure  sensors  arranged  in  a  500  m 
equilateral  triangle  (Figure  2,  bottom).  The  pressure  sensors  are  battery  powered,  internal 
recording  instruments  mounted  rigidly  near  the  sea  floor  inside  anchors.  The  autonomous 
instruments  were  synchronized  to  within  10  ms  using  accurate  temperature  compensated 
oscillators  and  thus  can  be  used  as  a  coherent  array  for  sea  and  swell  waves.  The  sample 
frequency  was  2  Hz.  Instruments  were  accurately  positioned  using  differential  GPS,  an 
underwater  acoustic  navigation  system,  and  direct  imderwater  tape  measurements.  Sensor 
spacing  uncertainties  range  from  about  0.4  m  for  the  shortest  array  legs  to  2.0  m  for  the 
longest  array  legs.  The  array  geometry  allowed  for  omni-directional  wave  measurements. 
Sensor  spacings  ranging  from  20  m  to  3(X)  m  were  chosen  to  accurately  resolve  swell  and 
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sea  waves  in  the  frequency  range  0.04  Hz  -  0.2  Hz.  Waves  with  frequencies  higher  than 
0.2  Hz  are  strongly  attenuated  at  the  sea  floor,  and  thus  not  accurately  measured  by  the 
pressure  sensors  (e.g.,  Herbers  and  Guza,  1991).  A  Datawell  Directional  Waverider  buoy 
was  deployed  within  the  array  perimeter  to  provide  measurements  of  the  high  frequency 
part  (0.2  Hz  -  0.5  Hz)  of  the  wave  spectrum.  This  surface-following  buoy  measures  3- 
component  accelerations  that  yield  wave  frequency  spectra  and  low-resolution  directional 
information  (e.g.,  Long,  1980).  Buoy  data  sampled  with  a  frequency  of  1.28  Hz  was 
transmitted  continuously  via  an  HF  radio  link  to  a  receiver  on  shore. 

The  nearshore  instrumentation  included  four  arrays  of  10  or  11  bottom-mounted 
pressure  sensors,  each  with  an  alongshore  aperture  of  210  m  and  a  cross-shore  aperture 
of  35  m  (Figme  2,  top).  The  arrays  were  positioned  at  offshore  distances  (relative  to 
F.R.F.  coordinates)  of  210,  260,  375,  and  500  m  in  nominal  depths  of  2.9,  3.7,  3.9,  and 
5. 1  m,  respectively.  The  nearshore  arrays  also  included  sonar  altimeters  that  provided 
accurate,  continuous  water  depth  measurements  (Gallagher  et  al.  1996).  All  instruments 
were  cabled  to  a  central  data  acquisition  center  onshore.  The  sample  rate  of  the  pressure 
sensors  was  2  Hz. 

Additional  directional  wave  data  were  available  from  a  permanent  array  of 
pressme  sensors  (maintained  by  the  U.S.  Army  Corps  of  Engineers  Field  Research 
Facility)  located  about  900  m  from  shore  in  8  m  depth  (Figure  2,  middle).  The  15- 
element  array  with  an  aperture  of  250  m  (alongshore)  x  120  m  (cross-shore)  and  a  2  Hz 
sample  rate  resolves  waves  with  frequencies  up  to  approximately  0.3  Hz. 

Bathymetry  at  the  site  is  characterized  by  a  gently  sloping  (*  1:250),  nearly  planar 
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inner  shelf  (Figure  1)  and  a  slightly  steeper  (« 1:100)  beach  (Figure  3).  Seaward  of  the 
20  m  array,  the  shelf  is  characterized  by  a  very  gentle  slope  (« 1:2500)  with  large  (scales 
0(5  km))  ridge-like  features  that  have  vertical  scales  of  0(5  m).  Refraction  computations 
show  that  the  effects  of  the  benign  shelf  topography  on  wave  transformation  across  the 
shelf  are  significant,  but  do  not  produce  strong  local  gradients  in  wave  energy 
(Hendrickson,  1996).  Daily  bathymetric  surveys  of  the  nearshore  region  (e.g..  Figure  3) 
indicate  that  the  beach  changes  over  the  course  of  the  experiment  were  small.  The  sand 
bar  located  300  m  offshore  (Figure  3)  remained  relatively  stationary  with  its  crest 
submerged  about  5  m  below  mean  sea  level,  and  did  not  strongly  affect  the  wave  shoaling 
process.  Close  to  shore,  the  morphology  was  more  dynamic,  with  the  development  of 
transient  bars  and  alongshore  depth  variations.  These  changes  occurred  well  shoreward 
of  the  nearshore  arrays,  and  thus  did  not  influence  this  study. 

In  the  following  analysis  of  wave  shoaling  evolution  from  20  to  2.5  m  depth, 
alongshore  depth  variations  are  neglected.  The  preliminary  analysis  presented  here  is 
focused  on  four  3-hour-long  case  studies  that  illustrate  the  observed  wave  shoaling 
evolution  in  both  non-breaking  and  breaking  conditions.  Prior  to  the  analysis,  the  data 
were  tested  for  stationarity  and  spatial  homogeneity  using  a  variety  of  methods.  Temporal 
changes  in  2-minute  average  wave  variances  were  used  to  ensure  stationarity  of  wave 
conditions  throughout  the  selected  3-hour-long  record.  Alongshore  homogeneity  over  the 
alongshore  extent  of  the  nearshore  arrays  was  verified  by  intercomparing  the  measured 
spectra  of  each  sensor  in  the  array,  as  well  as  the  coherence  and  phase  spectra  of 
redundant  alongshore  sensor  pairs.  For  example,  the  southernmost  sensors  of  the 
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nearshore  arrays  consistently  recorded  spectral  density  values  an  order  of  magnitude  lower 
than  other  array  sensors  for  waves  arriving  from  large,  oblique  southerly  directions. 
These  discrepancies  are  likely  caused  by  the  close  proximity  of  the  U.S.  Army  Corps  of 
Engineer’s  research  pier  (e.g.,  wave  blocking  by  pier  pilings),  and  therefore  data  from  the 
southernmost  sensors  was  not  used. 

Array  cross-spectra  were  calculated  based  on  Fourier  transforms  of  overlapping 
1024  s  segments.  Merging  1 1  frequency  bands  resulted  in  estimates  with  a  frequency 
resolution  of  0.0107  Hz  and  220  degrees  of  freedom.  For  each  frequency,  /,  the 
corresponding  directional  distribution  of  wave  energy,  S(0;/),  was  estimated  from  the 
array  cross-spectra  using  the  variational  technique  described  in  Herbers  and  Guza  (1990). 
The  direction  0  is  defined  relative  to  the  local  shoreline  orientation  with  0=0° 
corresponding  to  normal  incidence  to  the  beach  (waves  arriving  from  70"  true  N),  and  0 
is  positive  (negative)  for  waves  approaching  the  beach  from  northerly  (southerly) 
directions.  The  method  yields  the  smoothest  non-negative  distribution  function  that  is 
statistically  consistent  with  the  observations.  This  roughness  minimization  technique 
effectively  suppresses  spurious  peaks  in  S(0i/)  estimates  that  might  otherwise  be 
misinterpreted  as  multimodal  features  of  the  directional  wave  spectrum.  The  estimates 
of  S(0^  were  combined  with  an  estimate  of  the  surface  elevation  frequency  spectrum  E(/) 
(obtained  by  averaging  the  measured  auto-spectra  and  applying  a  linear  theory  depth 
correction)  to  form  the  wave  frequency-directional  spectrum  E(/,0)  (=  E(/)S(0i/)). 

Estimates  of  E(f,0)  were  initially  computed  at  all  array  locations  using  the  full 
array  geometry  (i.e.,  both  the  alongshore  and  cross-shore  legs).  Spurious  peaks  and  poor 
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convergence  of  S(0;/)  estimates  at  high  frequencies  observed  at  the  nearshore  and  8  m 
depth  locations,  suggest  that  nonlinear  and/or  wave  breaking  effects  contributed  significant 
errors.  The  estimation  technique  assumes  a  linear  dispersion  relation,  and  thus  is  sensitive 
to  nonlinear  changes  in  the  wavenumber  (e.g.,  amplitude  dispersion  and  bound  waves). 
Wave-breaking  induced  energy  losses  may  reduce  the  coherence  in  the  cross-shore  array 
legs  causing  additional  inversion  problems.  Frequency-directional  spectra  at  the  nearshore 
and  8  m  depth  locations  were  subsequently  re-computed  using  only  the  alongshore  legs 
of  the  arrays.  These  linear  array  estimates  are  insensitive  to  dispersion  errors  and  cross¬ 
shore  inhomogeneities,  but  cannot  distinguish  shoreward  propagating  waves  (-90‘’<9<90°) 
from  seaward  propagating  waves  (9O°<0<27O°).  Since  the  full  array  estimates  already 
confitmed  that  swell  reflections  from  shore  were  weak,  the  estimates  obtained  from 
the  linear  arrays  were  constrained  to  vanish  at  seaward  propagation  directions 
(9O°<0<27O’)  (Herbers  and  Guza,  1990).  The  resulting  estimates  were  nearly  identical 
to  the  full  array  estimates  in  the  energetic  part  of  the  spectrum  and  smooth  at  higher 
frequencies.  In  two  cases,  August  17  and  October  19,  estimates  of  E{40)  extracted  from 
the  20  m  depth  array  measurements  also  indicated  spurious  structure  and  convergence 
problems  at  high  frequencies.  The  estimates  were  recomputed  using  a  five  element,  linear 
sub-array.  In  the  August  17  case,  with  waves  arriving  from  southerly  angles,  the  northern 
leg  of  the  triangle  array  was  used  with  the  constraint  S(0;/)  =  0  at  angles  (-180°<9<-150% 
3O°<0<18O°).  Similarly,  for  the  October  19  case  with  waves  arriving  from  northerly 
directions,  the  southern  leg  of  the  triangle  array  was  used  with  the  constraint  S(9;/)  =  0 
at  angles  (-180°<9<-30’,  15O’<0<18O*). 
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Estimates  of  E(f,0)  in  20  m  depth  over  a  wider  frequency  range  (0.04  -  0.5  Hz) 
were  obtained  from  the  directional  wave  buoy  using  the  Maximum  Entropy  Method 
described  by  Lygre  and  Krogstad  (1986).  This  method  yields  a  miimodal  or  bimodal 
S(0iO  that  exactly  fits  the  measured  buoy  cross-spectra.  Directional  resolution  is  limited 
because  the  buoy  cross-spectra  define  only  four  moments  of  S(d]f).  The  buoy  spectral 
analysis  is  based  on  an  ensemble  average  of  54  data  segment  with  lengths  of  200  s.  The 
resulting  E(f,0)  estimates  with  0.005  Hz  frequency  resolution  were  smoothed  to  the  same 
resolution  as  the  array  estimates.  At  low  frequencies  where  they  overlap,  the  array  and 
buoy  estimates  are  in  good  agreement.  The  final  estimates  of  E(/j0)  in  20  m  depth  were 
formed  by  matching  the  array  estimates  in  the  frequency  range  0.04  -  0.2  Hz  with  the 
buoy  estimates  in  the  frequency  range  0.2  -  0.5  Hz. 
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m.  OBSERVATIONS 


Example  observations  of  the  shoaling  evolution  of  E(49)  are  shown  in  Figures  4 
-11.  The  four  cases  described  in  detail  below  span  a  wide  range  of  conditions,  including 
normal  and  large  obhque  mcidence  angles  and  both  non-brealdng  and  breaking  conditions. 
In  order  to  distinguish  linear  refraction  effects  from  nonlinear  harmonic  generation  and 
surf  zone  effects,  the  observed  evolution  of  E(f,9)  was  compared  to  simple  linear  theory 
predictions.  On  a  beach  with  straight  and  parallel  depth  contours,  the  transformation  of 
E(/j9)  is  governed  by 


C  C 

s  gs 

(Longuet-Higgins,  1957;  LeMehaute  and  Wang,  1982),  where  subscripts  d  and  s  denote 
wave  properties  at  a  deeper  and  shallower  cross-shore  location,  C  and  Cg  represent  the 
wave  phase  and  group  speed,  and  the  angles  9^  and  9^  are  given  by  Snell’s  Law 


C 

sin(9p=— sin(0  p 


(2) 


Eqs  (1)  and  (2)  were  used  to  transform  the  20  m  depth  estimates  of  E(f,9)  to  each  of  the 
shallower  array  locations.  The  resulting  predictions  are  included  in  Figures  4-11. 
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A. 


AUGUST  17 


The  August  17  conditions  are  characterized  by  a  steady  obhque-offshore  wind  with 
a  speed  of  8  m/s  and  direction  of  -140°.  Significant  wave  heights  were  less  than  0.5  m 
at  all  array  locations  and  breaking  occurred  in  the  region  shoreward  of  the  shallowest 
array.  Frequency-directional  wave  spectra  observed  at  the  six  arrays  are  compared  to  the 
linear  refraction  predictions  in  Figure  4.  Observed  and  predicted  directional  distributions 
of  energy  at  selected  frequencies  are  compared  in  Figure  5.  In  20  m  depth,  the  wave 
spectrum  in  the  dominant  0.15-0.30  Hz  range  is  broad  in  frequency  and  directionally 
narrow  with  waves  propagating  at  large  southerly  incidence  angles  nearly  parallel  to  the 
shore.  These  waves  are  strongly  refracted  over  the  inner  shelf  and  beach,  and  the 
observed  changes  in  wave  direction  are  well  described  by  the  linear  refraction  predictions. 
High  frequency  spectral  energy  levels  are  reduced  by  about  a  factor  5  between  the  deepest 
and  shallowest  sites,  consistent  with  predicted  spectral  levels  (left  panels).  The  observed 
refraction  of  the  smaller  and  broader  0.1  Hz  swell  peak  from  -50'  at  the  deepest  site  to 
-20”  at  the  shallowest  site  is  also  in  agreement  with  the  prediction  (Figiue  5,  left  panels). 
This  peak  is  amplified  during  shoaling  and  dominates  the  wave  spectrmn  at  the  shallower 
sites  (Figure  4),  in  agreement  with  the  prediction.  Whereas  observed  energy  levels  and 
mean  propagation  directions  are  in  good  agreement  with  the  predictions,  the  observed 
S(d/)  are  broader  than  the  predicted  extremely  narrow  S(Q  J)  (e.g.,  0.17  Hz  in  Figure  5). 
Possible  explanations  for  these  discrepancies  are  the  limited  directional  resolution  of  the 
array  and  alongshore  variations  in  wave  directions  caused  by  small  alongshore  depth 
variations  that  might  broaden  the  array  estimates  of  8(0;/). 
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B.  OCTOBER  2 


The  shoaling  evolution  of  a  more  energetic  wave  field  was  observed  on  October 
2  and  is  detailed  in  Figures  6  and  7.  Conditions  are  characterized  by  a  steady  onshore 
wind  with  a  speed  of  11  m/s  and  a  direction  of  55°.  In  20  m  depth,  E(f,9)  is  relatively 
narrow  in  both  frequency  and  direction  with  a  single  peak  at  approximately  /  =  0.17  Hz 
that  is  closely  aligned  with  the  wind  direction  (bottom  panels  of  Figures  6  and  7).  At 
higher  frequencies,  the  S(0;^  estimates  (extracted  from  the  directional  wave  buoy)  are 
bimodal  and  approximately  symmetric  about  the  wind  direction.  The  separation  of  the 
two  peaks  increases  with  frequency,  similar  to  recent  observations  by  Young  et  al.  (1995) 
and  Ewans  (1998).  Young  et  al.  (1995)  show  that  the  bimodal  stmcture  of  E(/j0)  at  high 
frequencies  can  be  explained  with  theoretical  predictions  of  transfers  of  energy  to  higher 
frequencies  through  nonlinear  wave-wave  interactions.  However,  the  8  m  depth  array 
estimate  of  E(f,0)  shows  high  frequency  waves  propagating  in  the  local  wind  direction 
and  suggest  that  the  bhnodal  distributions  estimated  from  the  buoy  data  may  be  erroneous. 
Errors  in  the  directional  buoy  estimates  of  S(0;/)  at  high  frequencies  may  be  introduced 
through  nonlinearities  in  the  buoy  response  or  the  inherent  lack  of  resolution  of 
directional  buoys  that  measure  only  four  integral  moments  of  S(Qif)  (e.g.,  Herbers  and 
Guza,  1990). 

The  observed  significant  wave  height  at  the  shallowest  array  is  about  10%  smaller 
than  the  predicted  value,  indicating  that  the  surf  zone  is  confined  to  a  region  shoreward 
of  the  shallowest  array  (upper  left  panel  in  Figure  6).  The  dominant  sea  peak  is  strongly 
refracted  over  the  inner  shelf  and  beach,  and  the  observed  changes  in  wave  direction  are 
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well  described  by  the  linear  refraction  predictions.  Similar  to  the  previous  case,  the 
predicted  refracted  peaks  are  directionally  narrower  than  the  observed  peaks  (left  panels 
in  Figure  7).  At  two  times  the  spectral  peak  frequency  (0.34  Hz)  the  E(/‘,0)  observed  at 
the  shallower  arrays  are  bimodal  with  one  peak  (0  ~  50°)  aligned  with  the  wind  direction 
and  the  expected  harmonic  peak  (0  «  20”)  aligned  with  the  primary  0.17  Hz  sea  peak 
(Figure  6).  This  harmonic  peak  is  strongly  amplified  during  shoaling  and  dominates  the 
high-frequency  spectrum  in  2.9  m  depth  (Figure  6). 

C.  NOVEMBER  1 

The  November  1  case  (Figures  8  and  9)  is  characterized  by  a  steady  oblique 
offshore  wind  with  a  speed  of  7  m/s  and  direction  of  -140”,  similar  to  the  August  17  case. 
Observed  significant  wave  heights  ranging  from  1.2  m  at  the  deepest  array  to  1.1  m  at 
the  shallowest  array  are  close  to  the  linear  theory  predictions  (left  panels  in  Figure  8) 
indicating  that  breaking  occurred  in  the  region  shoreward  of  the  shallowest  location.  In 
20  m  depth,  the  E{f,Q)  estimate  features  a  dominant  0.10  Hz  peak  with  a  large  southerly 
incidence  angle  (-60”).  The  high  frequency  part  of  the  spectrum  is  directionally  narrower, 
with  even  larger,  almost  shore-parallel  southerly  angles  (-80”),  similar  to  the  August  17 
case.  The  strong  refraction  over  the  inner  shelf  and  beach  of  the  dominant  wave  direction 
is  again  well  described  by  the  linear  refraction  predictions  (Figure  9,  left  panels).  As  in 
the  previous  cases,  predicted  directional  distributions  are  narrower  than  observed  at  the 
shallower  arrays.  At  high  frequencies,  second  (0.2  Hz)  and  third  (0.3  Hz)  harmonic  peaks 
are  observed  at  the  shallower  sites  which  are  approximately  directionally  aligned  with  the 
dominant  (0.1  Hz)  waves  (Figure  8).  Similar  to  the  October  2  case,  these  harmonic  peaks 
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are  amplified  between  5.5  and  2.8  m  depths  and  dominate  the  directional  distributions  at 
high  frequencies  at  the  shallowest  array  (Figure  9,  center  and  right  panels). 

D.  OCTOBER  19 

The  most  energetic  waves  in  the  SandyDuck  experiment  were  observed  on  October 
19  (Figures  10  and  11)  during  a  severe  Nor’easter  with  maximum  sustained  wind  speeds 
of  about  17  m/s.  The  surf  zone  extended  several  kilometers  offshore  with  a  significant 
wave  height  at  the  20  m  depth  array  of  about  3.9  m,  decreasing  steadily  to  1.56  m  at  the 
shallowest  array.  Propagation  directions  of  the  dominant  0.1  Hz  waves  are  nearly  normal 
to  the  beach  in  20  m  depth  and  virtually  unchanged  throughout  the  s\uf  zone  (Figme  1 1, 
left  panels)  even  though  spectral  energy  levels  are  reduced  by  more  than  an  order  of 
magnitude  at  the  shallowest  array  (Figure  10,  left  panels).  The  directional  distributions 
are  slightly  broader  than  predicted  with  linear  refraction  theory  at  the  shallower  arrays, 
but  the  broadening  of  8(9;/)  is  much  less  pronounced  than  that  observed  by  Herbers  et 
al.  (1998)  for  comparably  energetic  waves  breaking  on  a  shallow  bar.  At  twice  the 
spectral  peak  frequency  (0.19  Hz),  the  observed  spectra  show  a  transition  from  wind  seas 
in  20  m  depth  (propagating  in  the  local  wind  direction  0  «  80°)  to  second  harmonic 
waves  that  are  aligned  with  the  dominant  wave  direction  (9  «  0°).  In  8.2  m  depth  the  sea 
and  harmonic  peaks  are  approximately  equal  in  magnitude,  and  at  the  shallowest  arrays 
the  sea  waves  are  completely  submerged  in  the  harmonic  energy  levels.  A  similar 
transition  from  wind-generated  seas  to  harmonic  waves  is  observed  at  0.29  Hz  (Figure  11, 
right  panels).  As  in  the  August  17  case,  the  bimodal  8(0;^  observed  m  21.2  m  depth 
may  be  a  directional  buoy  artifact  rather  than  a  true  feature  of  the  wave  field.  In  8.2  m 
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depth,  a  refracting  sea  peak  from  40'  co-exists  with  a  developing  third  harmonic  peak  at 
5'  (aligned  with  the  dominant  waves).  This  harmonic  peak  broadens  during  shoaling  and 
dominates  the  spectrum  at  the  shallower  arrays. 
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IV.  SUMMARY 


The  evolution  of  the  frequency-directional  wave  spectrum,  E(f,0),  across  the  inner 
continental  shelf  and  beach  was  examined  with  extensive  measurements  collected  at  the 
Army  Corps  of  Engineer’s  Field  Research  Facility  near  Duck,  NC  during  the  recent 
SandyDuck  experiment.  The  barrier  island  field  site  is  characterized  by  a  straight 
coastline  with  a  gently  sloping  (1:250)  sandy  bottom  and  weak  alongshore  depth 
variations.  Arrays  of  bottom  pressure  sensors  were  deployed  on  the  inner  shelf  in  20  m 
depth  (5  km  from  shore)  and  on  the  beach  in  depths  ranging  from  2  -  5  m  (Figure  1). 
A  directional  wave  buoy  was  added  to  the  20  m  depth  array  to  measure  high  frequency 
seas  that  are  attenuated  at  the  sea  bed.  The  SandyDuck  arrays  were  complemented  by  a 
permanent  array  of  pressure  sensors  in  8  m  depth  maintained  by  the  Army  Corps  of 
Engineers.  High-quality  wave  data  were  collected  nearly  continuously  at  all  instrument 
locations  during  a  fom  month  period  (August  -  November  1997). 

The  preliminary  analysis  presented  here  is  focused  on  four  case  studies  that  span 
a  range  of  conditions  and  illustrate  the  observed  wave  shoaling  evolution  in  both  non¬ 
breaking  and  breaking  conditions.  Smooth  estimates  of  E(f,Q)  were  estimated  from  the 
cross-spectra  of  six  arrays  using  the  variational  technique  described  in  Herbers  and  Guza 
(1990).  In  20  m  depth,  the  array  estimates  of  E(^,0)  at  low  frequencies  were  combined 
with  buoy  estimates  at  high  frequencies  obtained  with  the  Maximum  Entropy  Method 
described  in  Lygre  and  Krogstad  (1986).  In  all  four  case  studies,  the  cross-shore 
evolution  of  the  dominant  wave  propagation  direction  agrees  well  with  simple  linear 
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theory  predictions  for  wave  refraction  on  a  beach  with  straight  and  parallel  depth 
contours.  In  the  more  energetic  cases,  the  observed  E(/',0)  show  the  development  of 
harmonic  peaks  at  the  shallower  sites  with  directions  that  are  approximately  aligned  with 
the  dominant  waves  (Figures  6,  8,  10),  qualitatively  consistent  with  theoretical  wave-wave 
interaction  rules  and  earlier  observational  studies  (Freilich  et  al.  1990;  Berbers  et  al. 
1998). 

In  both  non-breaking  and  breaking  conditions,  the  observed  E(/',0)  are  directionally 
broader  than  the  linear  theory  predictions,  in  particular  when  the  incident  wave  field  is 
directionally  narrow  (e.g.,  August  17,  Figures  4,  5).  These  discrepancies  may  reflect  the 
limited  resolution  of  the  arrays  and/or  alongshore  variations  in  wave  propagation 
directions  associated  with  weak  alongshore  depth  variations.  Whereas  Berbers  et  al. 
(1998)  observed  a  sudden,  dramatic  directional  broadening  of  wave  spectra  when  waves 
broke  on  a  shallow  submerged  sandbar,  the  present  observations  on  a  nearly  planar  beach 
do  not  indicate  a  strong  wave  breaking  effect  on  directional  wave  properties.  In  fact,  in 
the  most  energetic  case  (October  19)  when  the  surf  zone  extended  across  the  entire 
instrumented  transect,  the  directional  width  of  E(^,0)  remains  nearly  constant  upon 
shoaling.  These  differences  are  possibly  related  to  differences  in  nonlinear  shoaling 
evolution  and  wave  breaking  characteristics  (e.g.,  strongly  localized  breaking  on  a  sand 
bar  versus  gradual  breaking  on  a  planar  beach).  Further  work  is  needed  to  clarify  the  role 
of  the  sea  bed  profile  in  the  shoaling  evolution  of  directional  wave  properties. 
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APPENDIX 


Figure  1.  Plan  view  of  the  experiment  site  with  locations  of  the  nearshore,  8  m, 
and  20  m  arrays  and  associated  bathymetry  (depth  contours  are  relative  to  mean  sea 
level). 

Figure  2.  Detailed  plan  view  of  the  nearshore,  8  m,  and  20  m  arrays.  Pressure 
sensors  locations  are  represented  by  solid  dots  and  the  directional  wave  buoy  location  is 
indicated  with  a  circle. 

Figure  3.  Representative  beach  topography  during  SandyDuck  (comtesy  of  the 
Army  Corps  of  Engineer’s  Field  Research  Facility).  The  bathymetry  is  characterized  by 
a  small  shore-parallel  bar  located  275  m  offshore  and  nearly  straight  and  parallel  depth 
contours.  Bathymetry  changes  over  the  course  of  the  experiment  were  weak  with  the 
exception  of  the  dynamic  beach  face  region  shoreward  of  the  nearshore  arrays. 

Figure  4.  Wave  shoaling  evolution  from  19.8  m  depth  (bottom  panels)  to  2.5  m 
depth  (top  panels)  of  small  amplitude  waves  arriving  from  large  oblique  southerly  angles 
on  August  17,  1997.  Observed  (blue)  and  predicted  (red)  frequency  spectra  and 
significant  wave  heights  are  shown  in  the  left  panels.  The  color  contour  panels  show  the 
observed  and  predicted  frequency-directional  spectra  at  each  array  (units  cm^/Hz*;  0  =  0° 
corresponds  to  normal  incidence  to  the  beach)  with  the  water  depths  indicated  in  the 
upper  right  comer  . 

Figure  5.  Observed  (solid)  and  predicted  (dotted)  directional  distribution  of  wave 
energy  at  selected  frequencies  on  August  17,  1997.  Predictions  at  frequencies  higher  than 
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0.2  Hz  (only  measured  with  the  directional  wave  buoy  in  20  m  depth)  are  omitted. 

Figure  6.  Wave  shoaling  evolution  from  20.3  m  depth  (bottom  panels)  to  2.9  m 
depth  (top  panels)  of  moderately  energetic  waves  arriving  from  large  oblique  northerly 
angles  on  October  2,  1997  (same  format  as  Figure  4).  The  frequency-directional  values 
are  multiplied  by/^  to  enhance  the  directional  properties  of  high  frequency  waves. 

Figure  7.  Observed  (solid)  and  predicted  (dotted)  directional  distribution  of  wave 
energy  at  selected  frequencies  on  October  2,  1997  (same  format  as  Figure  5).  The  local 
wind  direction  is  indicated  with  a  dashed  vertical  line. 

Figure  8.  Wave  shoahng  evolution  from  20.8  m  depth  (bottom  panels)  to  2.8  m 
depth  (top  panels)  of  moderately  energetic  waves  arriving  from  oblique  southerly  angles 
on  November  1,  1997  (same  format  as  Figure  4).  The  frequency-directional  values  are 
multiplied  by/^  to  enhance  the  directional  properties  of  high  frequency  waves. 

Figure  9.  Observed  (solid)  and  predicted  (dotted)  directional  distribution  of  wave 
energy  at  selected  frequencies  on  November  1,  1997  (same  format  as  Figure  5). 

Figure  10.  Wave  shoaling  evolution  from  21.2  m  depth  (bottom  panels)  to  3.3 
m  depth  (top  panels)  of  large  amplitude  waves  arriving  from  directions  close  to  normal 
incidence  on  October  19,  1997  (same  format  as  Figure  4). 

Figure  11.  Observed  (solid)  and  predicted  (dotted)  directional  distribution  of 
wave  energy  at  selected  frequencies  on  October  19,  1997  (same  format  as  Figure  5).  The 
local  wind  direction  is  indicated  with  a  dashed  vertical  line. 
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